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ABSTRACT 

We re-examine scattering of photons near the Lya resonance in the intergalactic medium 
(IGM). We first derive a general integral solution for the radiation field around resonance 
when spin diffusivity is ignored. Our solution shows explicitly that recoil sources an absorp- 
tion feature, whose magnitude increases with the relative importance of recoil compared to 
Doppler broadening. This spectrum depends on the Lya line profile, but approximating it 
with the absorption profile appropriate to the Lorentzian wings of natural broadening accu- 
rately reproduces the results for a full Voigt profile so long as Tk ^ 1000 K in the IGM. This 
approximation allows us to obtain simple analytic formulae for the total scattering rate of Lya 
photons and the accompanying energy exchange rate. Our power series solutions converge 
rapidly for photons that redshift into the Lya resonance as well as for photons injected at line 
center. We confirm previous calculations showing that heating through this mechanism is quite 
slow and probably negligible compared to other sources. We then show that energy exchange 
during the scattering of higher-order Lyman-series photons can be much more important than 
naively predicted by recoil arguments. However, the resulting heating is still completely neg- 
ligible. 
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1 INTRODUCTION 

The radiative transfer of photons near the Lya resonance is cru- 
cial to understanding the high-redshift intergalactic medium (IGM), 
both because it determines the spin temperature of the 21 cm transi- 
tion l Wouthuvsen 1952; Field 1958) and because it affects the ther- 
mal history <Madau et alll997tlChen & Miralda-Escudel2004l) . 

The radiation field near this resonance has been examined 
a number of times in recent years. The earliest treatments ig- 
nored radiative transfer and assumed that the spectrum was fea- 
tureless around the line. Chen & Miralda-Escude 1 2004) were the 
first to solve (numerically) an approximate fo rm of the ra dia- 
tive transfer equation in this context (following lBaskdll98ll and 
iRvbicki & deirAntonidll994l) . They showed that, if photons red- 
shift toward the resonance, the spectrum develops an asymmetric 
absorption feature. As we will see explicitly below, the absorption 
feature is sourced by recoil in the scattering process: each scat- 
tering deposits an average energy AE = (hv a ) 2 /(m p c 2 ), where 
v a is the rest frequency of the Lya line. Thus photons lose en- 
ergy faster near the center of resonance, where they scatter more. 
To compensate for this increased "flow" speed, continuity requires 
that the amplitude of the background must decrease near resonance. 
This affects the scattering r ate of Lya photons and hence the spin 
temperature of the IGM. Hirata 1 2006) expanded on this method by 
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showing how to account for the hyperfine structure of the Lya line 
(see below). 

An alternative to the n umerica l approach of 

IChen & Miralda-Escudl 120041) and Iffiratal (2006) is to ap- 
proximate the spectrum analytically. This has a lon g history in 
resonant radiative transfer: Hummer &RvbickiliT99l summarize 
many of the advances. Of pa rticular interest to our problem is the 
treatment of lGrache\l ll989l) . who derived an analytic solution 
for the spectrum around a resonant transition when recoil is 
included. The analytic solution was obtained by approximating the 
absorption profile using the form appropriate for scattering in the 
Lorentzian wings provided by natural broadening. This assumption 
is valid when the optical depth is extremely l arge and the Doppler 
broade ning relatively small. Most recently, IChuzhov & Shapiro! 
(2005) rediscovered this solution and applied it to the problem of 
Lya transfer in the high-redshift IGM. In J3| we will show how 
these two kinds of solutions relate and evaluate when the analytic 
approximation is valid. We also compute the radiation field in such 
a way that the role of recoil becomes obvious. We examine the 
resulting total scattering rate in j3]and show that the approximate 
form proposed bv lChuzhov & Shapiro! i2005T) is a reasonably good 
match to the full numeric result. 

The line shape is also crucial for estimating the rate at which 
energy is transferred between the gas and the photon field. As de- 
scribed above, recoil during each scattering deposits some energy 
in the gas. If this were the sole mechanism for energy exchange, the 
IGM would rapidly be heated above the cosmic microwave back- 
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ground (CMB) temperature (Mada uetalll997l) . However, the ab- 
sorption feature actually cancels almost all of this heating. Consider 
a photon on the blue side of the line. This will be preferentially scat- 
tered by an atom moving away from the photon (so that it appears 
closer to resonance). The atom will then re-emit the photon isotrop- 
ically in its frame; in the IGM frame, the photon will therefore lose 
an energy ~ hAuo, where Avd is the Doppler width of the tran- 
sition. Photons that scatter on the red side, on the other hand, will 
tend to gain energy. The absorption feature develops so that this 
scattering "diffusivity" compensates for the recoil (i.e., so that more 
scattering occurs redward than blueward of the Lya transition). The 
net energy transfer is therefore much slower than naively expected 
JChen & Miralda-Escudel2004lRvbickl2006tlMeiksirl2006h . 

By employing their analytic approximation to the radiation 
field, IChuzhov & Shapiro! (2006 ) took a step toward finding a sim- 
ple solution for the net heating rate. In [J4] we take their approach 
further by deriving a fully analytic solution for heating by pho- 
tons redshifting into the Lya resonance as well as an approximate 
solution for photons injected at line center (either through recom- 
binations or cascades from higher hyn transitions). This allows us 
to examine how the heating rate varies with IGM temperature and 
optical depth. 

Of course, photons can redshift into any of the Lyn res- 
onances in the IGM. After a few scatterings, t hese p hotons 
are destroyed through cas cades to lower levels (Hirata 2006; 
IPritchard & Furla netto 200(f). The scattering rate is so small that 
recoil heating is negligible; however, all of the scatterings occur 
on the blue side of the line, so each depos its some fraction of 
the ato m's thermal energy in the gas as well. IChuzhov & Shapiro! 
1 2006) examined the analogous process in deuterium and found that 
it can be relatively strong. In Sj5| we show that the heating rate for 
Lyn photons is tiny even when frequency drift is included, because 
the photons scatter so far in the blue wing of the line. 

Hirata 1 2006) and Chuzhov & Shapiro 1 2005) examined how 
spin exchange affects the radiation spectrum. Because Lya 
transitions modify the ground- state hyperfine level populations 
jWoufhuvser]|l95l iFielj 1958), the photons can also increase or 
decrease their frequency during each scattering by an amount cor- 
responding to the energy defect of the 21 cm transition. This affects 
the flow rate of photons through the resonance (or more precisely 
the diffusivity) and hence the spectrum. However, because the level 
populations themselves depend on the Lya scattering rate, and be- 
cause the mean energy exchange per scattering depends on the level 
populations, including this effect on the spectrum requires an iter- 
ative solution. Fortunately, it is a small effect except at extremely 
low temperatures. Because we are most interested in the gross be- 
havior of the scattering and heating rates, we will neglect the spin 
diffusivity. When required, including it is relativ ely easy; th e steps 
are outlined in Furlanetto et al. (2006) (see also lHiratal20 06). 

In our numerical calculations, we assume a cosmology 
with Q m = 0.26, «a = 0.74, tt b = 0.044, and H = 
100ft km s _1 Mpc -1 (with h = 0.74), consistent with the most 
recent measurements l Sper gel et all2006l) . 



is (neglecting atomic recoil for the moment) 
1 dJ 



2 THE RADIATION FIELD NEAR THE Lya 
RESONANCE 

We let J be the comoving angle- averaged specific intensity (in units 
of photons per area per steradian). The equation of radiative transfer 



bj r 

4>(u) J + Hv a — + dv' R{v, v) J(u' 



in 



where nn is the hydrogen density, a a {v) = Xo^t") is the ab- 
sorption cross section, %a = (ive 2 /m e c)f a , f a is the absorption 
oscillator strength, and 4>{v) is the line profile. For our purposes, 
4> is given by the Voigt profile (which includes both collisional and 
natural broadening), 



a 



At 



a 2 + (z-i) 2 ' 



(2) 



with a = T /(An Avd), Y the inverse lifetime of the upper state, 
Avd/vq — (2ksTK /mc 2 ) 1 / 2 the Doppler parameter, Tk the gas 
temperature, and x = (v — vq)/ Avd the normalized frequency 
shift. The first term on the right-hand side of equation Q describes 
absorption, the second the Hubble flow, and the third re-emission 
following absorption. The redistribution function R(p, v') gives the 
probability that a photon absorbed at frequency v ' is re-e mitted 
at frequency v . The approximate form Ru(u, v') faenvevlll94lt 
lHummeJl 962). which assumes a Voigt profile with coherent scat- 
tering in the rest frame of the absorbing atom, is often used (see 
We must, however, also include recoil | Basko||l98 l|) and, for ex- 
act ca lculations, spin exchange l Hirata 2006; Chuzhov & Shapiro 
2005). The last term describes injection of new photons: C is the 
rate at which they are produced and tp(v) is their frequency distri- 
bution. 

This integro-differential equation simplifies considerably if 
we assume that the background spectrum is smooth on the scale 
of the average frequency change per scattering (which is Ax < 1; 
see Sj5)- If we neglect spin excha nge in this F okker-Planck approx- 



anck appr 
id 19941) 



imation, equation Q becomes iRvbicki & dell'Antoni 

^ {#"0^ + 2[ V 4>(x) + 7 ] J(*)} + Ci>{x) = 0, (3) 

wher e the Sobolev pa rameter 7 = Tq P , tgp is the total 
iGunn & Petersonl J 19651) optical depth of the Lya transition, and 
we have included at omic recoi l through the parameter 77 = 
(hv 2 )/(mc 2 Av D ) <Baskdll98ll) : this is the mean (normalized) 
frequency drift per scattering from recoil. Unfortunately, equa- 
tion |3j is not uniquely specified because there is some freedom 
in the drift and diffusivity imposed in the Fokker-Planck method. 
The form above matches that of IRvbicki & dell' Antonio! Il994l) ; 
however, it does not obey detailed balance, which requi res r\ —> 
V — l /(x + x a ) ~ n — l/x a where x a = v a /Avo iRvbicki 
2006). The correction is unimportant when IcbTk <JC hv a but is 
easily included in the analysis. Other form s of the Fok ker-Planck 
approximation have been examined by Meiksin (200(1). It is also 
straightforward to include the drift and diffusivity sourced by hy- 
perfine mixing in this formalism, so long as the spin temperature is 
known <Hiratal2006llChuzhov & Shapirol2005h . 

It is useful now to pause and note explicitly the scalings of the 
basic parameters of this problem; they will become useful later. We 
have Avd oc T^ 2 , so a oc T K 1/2 and T~ 1/2 . The Sobolev 
parameter 7 oc (1 + z)~ 3 ' 2 in the high-redshift limit. 

We will consider two sets of boundary conditions for equa- 
tion First, we let photons redshift into the resonance from large 
frequencies, with no injection term. To describe this we let Joo > 
be the specific intensity as x — > 00 and set C — 0. The second 
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Figure 1. Background radiation field near the Lyo resonance at z = 10, 
assuming a Voigt line profile. The upper and lower sets are for photons 
redshifting from infinity and photons injected at line center, respectively. 
(The former are normalized to Joo; the latter have J_oo = 1/2.) The solid 
and dashed curves take Tk = 10 and 1000 K, respectively. 

case allows injection at line center, so Cip(x) = CS(x), 1 and sets 
Joa = 0. In this case, we define J_oo to be the average inten- 
sity as x — > — oo as well. In either scenario, equation J3} is easy 
to integrate once, leaving us with a first order ordinary differential 
equation. 

The formal solution is most transparently obtained by chang- 
ing variables to iHummer & Rvbickil 1992b 



Ax' 



(4) 



+ 2(770 + 7)7 = 2^, 



(5) 



so that equation becomes 

dJ 
da 

where K = 7J00 for the continuous case, K — C for injected 
photons if x < 0, and K = for injected photons with x > 0. 
Obviously 



exp 



•2,1 



I 4>(a')da' + 2 7 cr 
Jo 



(6) 



is an integrating factor for this equation, from which the solution 
follows immediately. For injected photons with x > (so that K = 
0), it has the simple form 



J(x) = J(0) exp 



-2rjx — 2j 



Ax' 



(7) 



where J(0) is determined by continuity. 

A formal solution can also be written for K > 0, but in this 
case an alternate form is more physically illuminating. Here it is the 

1 Even if the initial Lyman-series absorption occurs well blueward of line 
center, the Lya photon that results from the cascade will be injected near 
line center because the atom passes through several intermediate states, each 
of which has a small natural width. 



Figure 2. Ratio of 8j in the "wing" approximation to the exact results (us- 
ing a Voigt profile). The solid and dashed curves assume a flat background 
spectrum and take Tk = 10 and 1000 K, respectively. The dotted and dot- 
dashed curves assume injection at line center, with Tk = 10 and 1000 K, 
respectively. In this case, we set 8j = J/J— oo for x > 0. All curves 
assume z = 10. 



absorption spike that is most interesting. To isolate its properties, 
we define Sj = (Joo — J)/ Joo', 2 note then that Sj > 0. The transfer 
equation takes the form 

dSj 



<t>-^; + 2{r 1 <t> + 1 )8j = 2770. 



(8) 



This has the same structure as the previous version, except that the 
sourcing term on the right-hand side depends on x. The same inte- 
grating factor yields the solution 



Sj(x) =2r) I dyexp 
'o 



-2777/ ~ 2 7 



da;' 



(9) 



This form makes it obvious that recoil sources the absorption spike. 
If the scattering was purely coherent, the gas and radiation field 
could not transfer any radiation and the spectrum would remain flat 
(see, e.g., Hummer & Rvbicki 1992). By sapping energy from each 
scattered photon, recoil increases the rate at which they redshift 
across the resonance. This increase in the "flow velocity" must be 
balanced by a corresponding decrease in the photon flux near the 
resonance. 

We show some example spectra in Figure ^ assuming that 
<j>(x) has a Voigt profile (see also lChen & M iralda-Escude 20041) . 
The upper curves assume that photons redshift into resonance from 
infinity; as expected, an absorption feature develops. It deepens at 
small temperatures, because, in that case, the energy lost from re- 
coil is large compared to the energy lost in each scattering (or 77 
is relatively large). The lower curves assume injection at line cen- 
ter. In this case, the spectrum spreads to large positive x when Tk 
decreases. 

2 For injected photons, Joo = 0, of course; then we make the substitu- 
tion Joo — > J_oo in the definition. We will see that Joo = J— 00 for a 
redshifting continuum. 
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Figure 3. Scattering integral as a function of IGM temperature. The thick 
solid, dashed, and dotted curves show (1 — S c ) for a Voigt profile at z = 
10, 20, and 30. The thin curves show the corresponding quantities using 
only the first-order term in equation 1171 . 



This nu merical solution is, of course , ide ntical to thos e 
presented by IChen & Miralda-Escudl J2004 and iHiratal <200ej) . 
once the appropriate line profiles, drifts, and diffusivities are in- 
ser ted. It is also a mor e gene ral form of the solutions provided 
by iHummer & Rvbickil ll992» (who neglected the recoil term) 
and|Chuzhov & Shapiro 1 2005). The latter made the approxima- 
tion (followingla^ that tj>(x) w 
a/ijvx 2 ), which is only accurate at \x\ ^> 1. We will refer to this 
as the "wing" approximation for convenience. This approximation 
allows the integrals over (ft -1 to b e performed analytically iGrachevi 
ll989tlChuzhov & Shapirol2005l) . For injected photons with x > 0, 
it is 



'GP 



The thick curves 



J(x) = J(0)exp -2rix 



2ty 72;" 
3 a 

while for a flat background or injected photons with x < 0, 



(10) 



5j(x) = 1r\ I dy exp 



3a vy 



3y 2 x + 3yx 2 ) 



2-qy 



.(11) 



This explains the discrepancy between the existing numeric and 
analytic results: the latter do not apply near the Doppler core of the 
profile. 

Figu re 13 shows t he ratio of the approxima te analytic form of 
iGrachevI <1989l) and lChuzhov & Shapiro! <2005l) to the exact spec- 
tra (computed with a Voigt profile). When the temperature is small 
(solid and dotted curves), the approximation is an excellent one. 
However, it begins to break down at large temperatures: for ex- 
ample, in the continuous case, it underpredicts 8j by ~ 10% at 
the center of the absorption spike when Tk = 1000 K. This is 
because the effective natural width decreases with temperature, 
so the thermal broadening becomes relatively more important in 
higher-temperature gas. For injected photons, the wing approxima- 
tion slightly shifts the curves to the right; when the decline at x > 
is sharp (as in the warm gas), the fractional deviation can be large. 
However, we find that the wing approximation is generally an ex- 



Figure 4. Scattering integral as a function of 7 

show (1 — S c ) computed numerically for a Voigt profile, while the thin 
curves show the corresponding quantities using only the first order term 
in equation 1171 . The dotted, short-dashed, long-dashed, and solid curves 
take Tk = 1, 10, 10 2 , and 10 3 K, respectively. The thin dot-dashed curve 
shows the approximate form proposed by Chuzhov & Shapiro 1 2005 ) for 
T K = 1 K. 



cellent one. In the next two sections, we will use this analytic form 
to study the scatterin g and heating rates, extending the approach of 
IChuzhov & Shapiro! <2005Ll200rt . 



3 THE Lya SCATTERING RATE 

The total rate at which Lya photons scatter (per hydrogen atom) is 

P a = ^ Xa / dvJ(y)(f>(u), (12) 



where J is now in proper units. Because each scattering can ex- 
change hyperfine states, this rate is crucial for determining the 
spin t e mperature of the 21 cm transition in the IG M ( Wouthuv^er 
|1952 }: iFiel d [ 1958 ]: iMadau et alJ fl997^ IChen & Miralda-Escudd 



boO^lHiratalboOcf IChuzhov & Shap!rll2005l) . T he Wouthuysen" 
Field coupling strength can be written as (e.g., iFurlanetto et alJ 
2006 ) :i 



(13) 



_ 167TXnJoo T* „ 

27A 10 T 7 



where A10 = 2.85 x 10~ 15 s" 1 is the spontaneous emission coef- 
ficient of the 21 cm transition, T* = 0.068 K is the energy defect 
of that transition, T 7 is the CMB temperature, and 



dx 4>{x) 



(14) 



depends only on the shape of the background spectrum. Note that 
S a < 1, because recoil always induces an absorption feature. 



For injected photons, one must substitute Jo. 
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In general, S a must be computed numerically; even in the 
wing approximation, there is no closed-form analytic solution. 
However, recall that 4>{x) is sharply peaked around x = 0, while J 
varies slowly near resonance (even in the injected case). Thus we 
can approximate J w J(0) everywhere inside the integral; from 
the normalization of d> we thus have 



S a » Sj(0). 



(15) 



In the wing approximation, this is easily computed from equa- 
tion mv 



1-S a 



4a 



+(3^)^(1; 1^; 

' 3 3 27 



(16) 



4?i 



8tt 



-a 2 + -a 3 + ...,(17) 
3^3r(2/3) 3>/3r(l/3) 3 

where Bi(a;) is an Airy function, 1F2 is a hypergeometric function, 
and 



a — r) 



( 3a 
I 2tt7 



1/3 



0.717T, 



-2/3 



10" 



1/3 



(18) 



where Tk is in degrees Kelvin. When a is small, we therefore have 
(1 — S a ) oc T^ 2 ^ 3 Tqp . This scaling gives some intuition for how 
the coupling strength varies in the IGM. As in Figure Q the ab- 
sorption spike becomes less and less significant as Tk increases; 
thus we must have S a — » 1 (its value without recoil) in a warm 
IGM. The perturbation increases with optical depth because that 
increases the number of scatterings (and hence the energy loss due 
to recoil). 

We show the dependence of (1 — S a ) on temperature in Fig- 
ure |3] and the dependence on the Sobolev parameter (or optical 
depth) in Figure|4| The thick curves show the numeric solution for a 
Voigt line profile and for a continuous background spectrum, which 
we denote S c . The case with photons injected at the line center has 
a nearly identical scattering integral, because Sj(0) is identical in 
the two cases; only at high temperatures does the structure around 
resonance matter. The thin curves show the first-order (in a) ap- 
proximation of equation i 171 . We see that this provides an excellent 
match at Tk > 10 K, especially when 7 is relatively large (i.e., at 
lower redshifts). 

Note that we have actually made three approximations here: 
(i) a constant J across the line; (ii) the wing approximation; and 
(iii) the small a approximation. The culprit at small Tk is the third. 
Here a is large and the power series approximation breaks down. 
However, even including just terms up to a 3 dramatically improves 
the estimate, with errors < 10% so long as Tk > 2 K. This demon- 
strates that the first approximation is an excellent one here: at such 
small temperatu res, 4>(x) is extremel y sharply peaked. The second 
is equally good. Chuzhov & Shapiro (2005) proposed the fit 



8j(0) « 1 - exp(-1.12a), 



(19) 



which retains the first order behavior of Sj (0) at small a (and hence 
is reasonably accurate) and fits the behavior for a ~ 1 much bet- 
ter. The thin dot-dashed curve in Figure|4]shows how well this ap- 
proximation does at Tk = 1 K; it typically differs from the exact 
solution by ~ 5%. 

The overall agreement worsens at large temperatures as well. 
Here a is small, so the power series in equation 1 1 71 converges 
rapidly and approximation (iii) is excellent. The problem lies in- 
stead with the other two. As we have seen, the wing approximation 



breaks down once Tk exceeds ~ 1000 K. This causes up to a 10% 
underestimate of <5j(0). At the same time, approximation (i) breaks 
down and the region around resonance starts to contribute to the 
scattering integral. This causes a < 20% overestimate of (1 — S c ) 
compared to the exact result; fortunately, these two effects partially 

cancel. 

In summary, the fit proposed by Chuzhov & Shapiro 
(in eq.[T9l is an excellent approximation to S a unless high accuracy 
is required. However, we emphasize that, in order to inclu de spin 
transfer properly, one must still use an iterative procedure iHiratd 
120061 IChuzhov & Shapirol2005t) . 



4 HEAT EXCHANGE FROM Lya SCATTERING 
4.1 Continuous Background 

The rate at which the radiation field deposits energy in the gas (per 
unit volume) is <Chen & Miralda-Escudel2004) 



r 



dv (Joo — J), 



(20) 



where we have assumed 1/ w v a across the absorption feature and 
J is again in proper units. The physical interpretation of this form is 
straightforward: in the absence of scattering, the absorption feature 
would redshift away to infinity. To keep it in place, the photons 
must lose energy at the rate given in equation 1201 . More formally, 
it can b e derived from the average energy exchange per scattering 
IChuzhov & Shapirol2 006 ) through integration by parts and the use 
of equation |3j- 

Thus the heating rate depends on 



h = 



= 277 



dx 5j(x), 
dye" 2 ™ 



da; exp 



dx' 



(21) 



(22) 



This cannot be done in closed form for an arbitrary line profile, but 
the accuracy of the wing approximation makes it extremely useful 
in understanding the solution. 4 In this case, both integrals can be 
done analytically, yielding 



(!) 



-1/6 



la = 

where 



7 ' T7 J 



.3/2 



1/3 



/?[Ai 2 (-/3) + Bi 2 (-/3)] , (23) 



0.99T, 



-2/3 



K 



7 

10 6 



1/3 



(24) 



and Ai(x) and Bi(x) are the Airy functions and Tk is in degrees 
Kelvin. Note that this solution is exact within the wing approxima- 
tion. 

We can again find a simple and useful approximation by ex- 
panding in powers of fi. We find 

<2^\ 5 / 3 (a^ 1/3 r 
. 1 , 



3>/ , (f) . 

3 2/3 /3 3 
+ r 2 (l/3) + ' 







r 2 (2/3) 



r(i/3)r(2/3) 

(25) 



4 Note that Chuzhov & ShaDirJ l2O06tl apparently calculated the heating 
rate numerically in the wing approximation. 
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Figure 5. Heating integral for continuous injection. Left panel: The thick solid, dashed, and dotted curves show I c for a Voigt profile at z = 10, 20, and 30. 
The thin curves show the corresponding quantities using only the first order term in equation 1251 . Right panel: Ratio of the power series approximation to I c 
(using the wing approximation) to the exact value. The thick and thin curves retain terms to order j3 and /3 3 , respectively. 



Thus, we see I c oc T A , 5/6 7 _2/3 ; because 7 = r^p oc (1 + z)~ 3/2 
at high redshifts, we expect I c oc (1 + z) at fixed temperature. The 
heat input per atom per Hubble time (at constant Joo) is therefore 
AT oc HAv D I c /(n m H) oc T~ 1/3 (l + z)~ 2 . These scalings are 
close to those estimated by Chu zhov & Shapirol fcOOfJ) from their 
numerical results. 

We show our solution for I c as a function of Tk in the left 
panel of Figure[5]and as a function of 7 in Figure|6| In each of these 
panels, the thick curves use the full Voigt profile, while the thin 
curves use the first-order term (in /3) of equation I25i : we expect 
the latter to be valid when Tk > 10 K. The right panel of Figure|5| 
shows the ratio of the approximate and exact solutions; here the 
thick curves retain only the lowest order term, while the thin curves 
include terms up to A : these are necessary for Tk < 10 K. As be- 
fore, the expansion in equation 1251 converges rapidly at higher 
temperatures, but the wing approximation begins to break down. 

Obviously the predicted scalings are reasonably accurate; 
the heating rate decreases with temperature (because recoil is 
relatively inefficient) and increases with tgp (along with the 
scattering rate). The higher-order terms, and the Voigt profile, 
slightl y decrease the d ependence on these parameters. As shown 
bv lChen & Miralda-EscudJ 120041) . Lya heating is probably slow 
compared to other processes, and the wing approximation (in the 
full analytic expression for small temperatures and the power se- 
ries form otherwise) should be adequate for most purposes. 



100 - 




Figure 6. Heating integral for continuous injection. The thick curves 
show I c computed numerically for a Voigt profile, while the thin curves 
show the corresponding quantities using only the first order term in equa- 
tion 1251 . The dotted, short-dashed, long-dashed, and solid curves take 
T K = 1, 10, 10 2 , and 10 3 K, respectively. 



4.2 Injection at Line Center 

For photons injected at the li ne center, a similar exercis e shows that Again, we work in the wing approximation to gain some intuition, 

the relevant integral is IChen & Miralda-Escudel20o4) The second integral can be written in closed form; the first is 

fQ f oc t / \ />0 I /* oc , , . 

J ( x > ,™ / 1 a d v ""t. .3 



«!«*,(«)-/ d*^- (26) / d»fc = -j=J-J o ^^{^-^y) 
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Figure 7. As Fig.|5] but for injection at line center. In the right panel, the approximate versions include terms up to order p (thick curves) and /3 2 (thin curves). 



xerfc 



7T7 

la 



(27) 



Unfortunately, the complementary error function prevents a closed 
form solution. However, note that the exponential term implies that 
the integral is dominated by the region where the argument of the 
error function is small. Expanding it to lowest order, we then obtain 
a power series solution in (3: 



li 



1/3 cc 

E 

i=0 



(28) 



The first few terms have {Ao, Ai, A2) = 
(-0.6979, 2.5424, -2.5645). Retaining only the zeroth-order 
term, t he scaling wit h 7 an d Tk is again close to that proposed by 
IChuzhov & ShapirdfcOOfj) at T K > 100 K. 

Figures and [S| show It for the same parameters as in Fig- 
ures [5] and [6] again we compare the approximate form with the 
exact solution (including the full Voigt profile). In this case the 
dependence on both Tk and 7 is considerably more complicated. 
Most interestingly, injected photons can both heat the gas (when 
Tk < 10 K) and cool it. Physically, cooling can occur because 
more photons scatter on the red than the blue side of the line; in 
such events, the re-emitted photon generally has a higher energy 
in the IGM frame and so removes heat from the gas. In the high- 
temperature regime (Tk ^ 100 K), the cooling rate falls slightly 
when 7 decreases and when Tk increases. At small temperatures, 
the exchange switches to heating because the feature is so broad 
compared to the Ax ~ 1 frequency change per scattering. 

We also show the approximate form (eq. !28> in these panels 
(note that we must include 0° and /3 1 terms). It is substantially less 
accurate at first-order in j3, only approaching the exact solution at 
Tk > 200 K; shortly thereafter, the Voigt profile becomes signifi- 
cant. However, the thin curves in the right panel of Figure0show 
that carrying the series expansion to /3 2 is quite accurate throughout 
the range Tk ^ 10 K. In the injected case, the wing approximation 
is less useful because no analytic solution exits. Thus we recom- 



„- 20 




Figure 8. As Fig.|6| but for injection at line center. Here the approximate 
version include terms up to order j3. 



mend numerical integration of equation I26i when high accuracy is 
required (especially at small temperatures). 



5 SCATTERING OF LY7Y PHOTONS 

Consider a photon that redshifts into a Lyn line (with frequency 
v n ) at redshift z r \ its frequency at redshift 2 is therefore i/ z = 
^ n [(l + z)/(l+2 r )]. The accumulated optical depth it has traversed 
by that point is 
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-(*) = 



dz 



a(u z )nm(z)c 



(29) 



(l + z)H(z) ' 

where a refers to the line of interest. The lGunn & Petersonl 1 19651) 
optical depth is of course the total optical depth experienced by 
such a photon, orr(z <C z r ). 

We are interested in determining the surface at which such a 
photon will first scatter. For the extremely optically thick Lyman- 
series lines of the IGM, this first scattering occurs far in the wings 
of the line, so we can set <p(x) « a/(irx 2 ). Further assuming the 
high-redshift limit [H(z) oc (l+z) 3/2 ] and letting Az = z-z r < 
z r , we have 



-(*)■■ 



a (1 +z, 



HoV^o ^ &z v 2 



S5/2 



(30) 



where Xn is evaluated for the line of interest. Re-expressing Az in 
terms of a;, we find 



x(t) 



1650 



tT, 



1/2 



m' fe) m ( 



'l + Zr 

20 



3/2 



(3D 



where we have normalized v, the spontaneous emission coefficients 
A n , and the inverse lifetimes F n to the values appropriate for Ly/3. 
Higher Lyman-series photons have significantly longer lifetimes 
and hence scatter nearer line center; for example, Lye photons have 
a coefficient ~ 18. 

By setting r = 1 in equation 13 11 . we see immediately that the 
first scattering occurs well blueward of resonance; we will denote 
this location x 1 . Lya photons cannot be destroyed during scattering 
(except, of course, in the exceedingly unlikely event that a collision 
occurs while the atom is excited), so this first scattering limit has lit- 
tle physical interest. However, higher Lyman-series photons can be 
destroyed, because the excited state can cascade to an intermediate 
level. The de struc tion probabilities per scattering are compiled by 
Iffiratal <200d) and Pritchard & Furlanetto 1 2006); they are ~ 10% 
for hy/3 and ~ 20% for higher-level transitions. Thus, each such 
photon scatters only a few times before vanishing; so long as they 
remain in the wings, the fcth scattering will occur at Xk ~ xi/k. 

Because the photons are far out on the blue wing during each 
of these scattering events, they will deposit some fraction of their 
energy in the gas, heating it slightly. Our next goal is to calculate 
the net energy exchange with the IGM as these photons scatter 
and eventually disappear. We begin with the redistribution func- 
tion Rn(x, x'), which gives the probability that a photon absorbed 
at frequency x' is re-emitted at frequency x, assu ming co herent 
scattering in the rest frame of the absorbing atom {Henvevl lT94lL 
lHummeJl962l) . thus ignoring recoil: 



R n(x,x) = ^ 



- tan" 



due u 

x—x\/2 
1 (X — U 1 

a 



tan 



(32) 



where x = max(x, x') and x = min(a;, x'). This is normalized so 
that 



da; Rn(x, x). 



(33) 



In our case (x' ^> 0), the redistribution function is sharply 
peaked around x . We thus write x — x' + Ax and expand the 
inverse tangents to third order in (u, Ax) about x (note that, al- 
though u can be arbitrarily large, the exponential guarantees that 
only small values contribute to the integral). Equation <32> then be- 



comes 

Ru(x,x') 



T 3/2 X ,2 



due' u [(2m -Ax) 



+ 



(Aa 



Ax\/2 

2uAx 



X > X 



(34) 



with a similar expression when x < x'. We are interested in the 
mean energy loss in each scattering, 



/oo 
da; Aa; Rn(x, x') 
- oo 



-2 

rV2r 



dAa; Aa; 



due u (2m — Aa 



(35) 



(36) 



Ax/2 



where we have used the normalization of Ru and substituted our 
expansion for the redistribution function in the second equality. The 
integrals are elementary and are most easily performed by switch- 
ing the order of integration; the simple result is 

(Ax\x) = -l/x'. (37) 

This is identical to th e second term in equation (3) of 
IChuzhov & Shapiro! 12005T) in the appropriate limit. 

As expected, photons tend to lose energy to the gas, but only 
slowly. Physically, we have assumed that the scattering is coher- 
ent in the rest frame of the atom. Thus, in the IGM frame, the gas 
tends to gain energy if the scattering atom travels away from the ini- 
tial photon and to lose energy if the atom travels toward it. When 
scattering occurs blueward of resonance, the former have a slightly 
higher cross section for absorption because of the small blueshift 
imparted to them by their thermal velocity; thus the net effect is en- 
ergy transfer to the gas. However, far out on the wings of the line, 
the difference in cross sections from this displacement is small, and 
the heating is weak. 

In contrast, consider scattering with zero natural width (or in 
other words where the total optical depth is small, so that the initial 
scattering occurs in the Doppler core). In that case, Ru simplifies 
to Ri(x,x') = erfc(|x|)/2 (again assuming isotropic sc attering in 
the rest frame), where \x\ = max(|a;|, |x'|) (Unno 1952). This has 
a flat core between (—x r , x'), so the typical energy lost in the initial 
scattering is ~ x . In this case, photons on the blue side are only 
scattered by atoms moving away from them (so that their frequency 
lines up with the resonance), and in the lab frame the re-emitted 
photon typically shifts by a full Doppler width. Thus, in the few- 
scatterings limit, heating will be most efficient inside the lin e core. 
This is the case considered (for deuterium) by Chuzh ov & Shapirol 
(2006). 

Returning to the Lyn lines, where all the interactions occur in 
the wings, the net frequency shift in s scattering events is 

s(l + s) 



Ax to 



2xi ' 
O.OSST^ 2 



l(l + S r , 



110 

l + Zr \ ~ 3 ' 2 

20 



Vn\ ( A@ 



(38) 



(39) 



Again, we have normalized to the values appropriate for Ly/3 pho- 
tons. Of course, we have assumed that the scatterings occur in the 
wings of the line. If the accumulated drift carries the photon to- 
ward line center, subsequent scattering will occur symmetrically 
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and heating will be negligible. Thus we must have Aa;tot < x\. 
This is marginally true for Lye, which has s e = 5 and a coefficient 
0.77T^/ 2 when appropriate values are inserted into equation 1391 . 

It is useful to compare this drift to that due to recoil itself, 
which we ignored by using Ru. This has As rcco ii = r\ per scatter- 
ing. Thus 




the coefficient is 2.35 for Lye photons. Thus, at reasonably large 
temperatures, the frequency drift from repeated scattering over- 
whelms recoil. However, unless the IGM is already warm , it never 
dominates by a large factor. IPritc hard & Fu rlanetta 12006) showed 
that recoil provides a negligibly small contribution at all tempera- 
tures. Thus, in practice heating from Lyn scattering is never signif- 
icant: only if Tk 2> T 7 could it possibly matter, but in that case 
other, much stronger, heating agents must already be present. 

As we have shown, Aa; is largest (~ 1) w hen scatte ring occurs 
near th e line core. One example, considered by Chuzhov & Shapiro 
( 2006), is the deuterium Ly/3 resonance, for which the optical depth 
is of order unity. Then the energy transfer is much more efficient. 
This situation has different temperature dependence. In either case, 
AT oc hAvo (Ax) per scattering. When absorption is in the line 
center (Ax) ~ 1; in the wings, we have seen that natural broad- 
ening controls the cross section and (Ax) oc 1/xi oc Aun- So 
AT core oc TjJ/ 2 and AT w i ng oc Tk. Deuterium Ly/3 turns out to be 
the most important transition (aside from hydrogen Lya) in heat ex- 
change. Of course, this energy is injected into the deuterium, rather 
than the hydr ogen, to which it must be transferred by collisions. 
According to Chuzhov & Shanird 120051) . this is relatively ineffi- 
cient, so Lya heating still dominates by a large factor. As a result, 
the deuterium temperature may become quite large (T ~ 10 4 K), 
where the wing approximation breaks down and the full Voigt pro- 
file must be used. 



6 DISCUSSION 

We have examined both analytic and numeric solutions for the ra- 
diation field near the Lya resonance and used them to compute 
the total scattering rate and the IGM heating (or c ooling) rate. We 
sho wed that the approximate a nalytic solution of iGrachevl i 19891) 
and lChuzhov & Shapiro! 120051) . in which scattering in the wings 
dominates, is accurate so long as Tk < 1000 K. At higher tempera- 
tures, thermal broadening becomes important. Fortunately, the scat- 
tering correction S a ~* 1 at large temperatures. So the approximate 
fit presented by Chuzhov & Shapiro 1 2005) - our equation I19i - 
turns out to be reasonably accurate (to several percent) whenever 
Tk ^1K, For higher accuracy, equation < Kit can be used. 

We then used this analytic solution to examine the heating (or 
cooling) from the scattering near line center. For the case of pho- 
tons that redshift toward the resonance, we obtained a fully ana- 
lytic solution (in terms of Airy functions) under the approximation 
that all scattering occurs in the wings. The arguments of the Airy 
functions are typically small, so a power series expansion is illumi- 
nating; it shows that the heating rate per atom and per Hubble time 
is proportional to T^ 1/,3 (l + z) -2 ; this is an excellent approxima- 
tion at Tk > 10 K. In the case of photons injected at line center, we 
obtained a power series solution that converges reasonably rapidly. 



The lowest order term is reasonably accurate for Tk > 100 K, and 
in this regime the cooling rate per atom and per Hubble time is pro- 
portional to Tj/ 3 (1 + z) _5//2 . Photons injected in this way only 
heat the gas when T K <WK. 

We have neglected drift and diffusivity caused by spin ex- 
change during Lya scattering. This is relatively e asy to inco rporate 
into the approximate solution IChuzh ov & Shapiro 2005). How- 
ever, the spin temperature actually depends on the spectrum around 
Lya, so an i terative process must be e mployed when Tk is sm all 
lHiratal2006l) . The steps are outlined in lFurlanetto et alj f2006). 

As a final thought, it is useful to estimate the heating rate from 
Lya scattering to gauge its importance relative to other processes. 
For simplicity we will consider continuous injection (i.e., photons 
redshifting into the Lya resonance). Inserting our lowest order ap- 
proximation for I c (from eg. 125 \ into equation 1201 . we find 

2 e a 0.80 scg / JO_ \ 

3 Hnmk B T K ~ r 4/3 S a \l+z)' 

where the left hand side is the fractional temperature change per 
Hubble time. On the right hand side, we have rewritten J{y) in 
terms of the 21 cm coupling efficiency x a (see eq. H3> . The 21 cm 
spin temperature departs from the CMB temperature when x a ~ 1, 
so this is a convenient gauge for the background fluxes at which 
heating is relevant. Clearly, Lya heating is negligible unless the 
initial temperature is also small. (Note that this approximation for 
7 C overestimates the heating at low temperatures, so the actual heat- 
ing is even smaller than predicted by ea.1411) Because, e ven without 
any heating, T K = 2.5 K at z = 10 ISeaeer et alJl999h . Lya scat- 
tering is unlikely to be significant in this context. Other proces ses, 
especially X-ray heating, are probably much more important lOhl 
1200 it iGlover & Bran j20o3: lFurlanettol2006l) . 

We thank G. Rybicki and M. Furlanetto for helpful discus- 
sions. 
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